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ABSTRACT 

Highly-ionized fast accretion-disk winds have been suggested as an explanation for a variety 
of observed absorption and emission features in the X-ray spectra of Active Galactic Nu- 
clei. Simple estimates have suggested that these flows may be massive enough to carry away 
a significant fraction of the accretion energy and could be involved in creating the link be- 
tween supermassive black holes and their host galaxies. However, testing these hypotheses, 
and quantifying the outflow signatures, requires high-quality theoretical spectra for compar- 
ison with observations. Here we describe extensions of our Monte Carlo radiative transfer 
code that allow us to generate realistic theoretical spectra for a much wider variety of disk 
wind models than possible in our previous work. In particular, we have expanded the range 
of atomic physics simulated by the code so that L- and M-shell ions can now be included. 
We have also substantially improved our treatment of both ionization and radiative heating 
such that we are now able to compute spectra for outflows containing far more diverse plasma 
conditions. We present example calculations that illustrate the variety of spectral features pre- 
dicted by parameterized outflow models and demonstrate their applicability to the interpreta- 
tion of data by comparison with observations of the bright quasar PG121 1+143. We find that 
the major features in the observed 2-10 keV spectrum of this object can be well-reproduced 
by our spectra, confirming that it likely hosts a massive outflow. 

Key words: radiative transfer - methods: numerical - galaxies: active - X-rays: galaxies - 
galaxies: individual: PG121 1+143 



1 INTRODUCTION 

The vast energy released in accretion by a supermassive black hole 
means that active galactic nuclei (AGN) could have a fundamental 
role in the formation and evolution of galaxies (see Cattaneo et al. 
2009). This depends however on the nature of the physical mech- 
anisms that couple the AGN to its host galaxy. One such possible 
AGN feedback mechanism is that of a massive outflow, launched 
somewhere close to the accreting black hole. It has already been 
shown that such flows might be able to account for observed cor- 
relations between properties of nuclear black holes and their host 
galaxies (see e.g. King 2003, 2005). 

Direct evidence for AGN outflows and quantification of their 
properties, however, requires interpretation of observational data. 
Perhaps the best evidence for energetically important outflows 
around AGN comes from the detection of X-ray absorption lines 
identified with significantly blueshifted line transitions of highly 
ionized material (see Turner & Miller 2009 for a recent review of 



X-ray observations of AGN). The most promising origin for mass- 
loss in the immediate vicinity of an accreting supermassive black 
hole is an accretion disk wind. Hydrodynamical simulations have 
demonstrated that such flows are plausible (e.g. Proga & Kallman 
2004). However, these simulations also suggest that disk winds 
are likely to be very complex and thus to interpret their spectro- 
scopic signatures quantitatively requires detailed synthetic spectra 
for comparison. 

In previous work (Sim 2005, Sim et al. 2008, hereafter 
Paper I), we developed a numerical code for computing syn- 
thetic X-ray spectra for outflow models. We showed that simply- 
parameterized disk wind models could readily account for strong 
blueshifted absorption features associated with Fe XXV and 
Fe XXVI and we explored some of the effects of the outflow den- 
sity, geometry and ionization state on these absorption line proper- 
ties. Our models also confirmed that a very highly ionized outflow 
could affect the X-ray spectrum in more subtle ways. In particular, 
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the blueshifted Fe absorption lines have associated emission fea- 
tures that are broad and can develop extended red-skewed wings 
owing to the effects of electron scattering in the flow (see also 
Laming & Titarchuk 2004; Laurent & Titarchuk 2007). We con- 
cluded that, although complex, outflow models have the promise 
to simultaneously account for several of the observed X-ray spec- 
troscopic features of AGN. We illustrated this via a direct com- 
parison to observations of the well-known narrow line Seyfert 1 
galaxy Mrk 766. However, that work was limited to the treatment 
of only the atomic physics necessary for the most highly ionized 
atomic species (specifically, He- and H-like ions). Although AGN 
provide copious ionizing radiation, the region responsible for the 
primary X-ray emission is assumed to be small and centrally con- 
centrated. Therefore more distant portions of an outflow - or those 
shielded from the X-ray source by other material - need not be 
so significantly ionized. Thus a more realistic study requires that a 
wider range of ionization conditions can be considered. In partic- 
ular, photoelectric absorption and/or fluorescent emission by lower 
ionization states of Fe may affect the X-ray spectrum, both in the 
Fe K region and at lower X-ray energies. For lower typical ioniza- 
tion conditions, the relative importance of the lighter elements also 
grows significantly and they can lead to strong photoelectric ab- 
sorption for energies ~ 1 keV and associated discreet line features 
in the soft X-ray spectrum as have been reported in observations 
(e.g. Pounds et al. 2003, 2005; Pounds & Reeves 2007, 2009). 

Here, we extend our Monte Carlo radiative transfer code de- 
scribed in Paper I to treat the physics of L- and M-shell ions of 
astrophysically abundant elements. We also substantially improve 
the treatment of the ionization state and introduce a simple means 
of estimating the kinetic temperature in the outflow, thereby elimi- 
nating this as a free parameter of the models. In addition, we modify 
the means by which the orientation-dependent spectra are extracted 
from our simulations in order to suppress the level of Monte Carlo 
noise. These improvements to the code and atomic data used are de- 
scribed in Section 2. In Section 3 we demonstrate that our improved 
treatment of ionization leads to quantitatively good agreement with 
well-known ID codes and in Section 4 we discuss results from a 
calculation made with the improved code. Finally, in Section 5, we 
discuss the comparison of spectra computed with our models to ob- 
servations of the quasar PG121 1+143 before drawing conclusions 
in Section 6. 



2 METHOD 

2.1 Code overview 

The code operates by performing a sequence of Monte Carlo radia- 
tive transfer simulations employing an indivisible packet scheme 
(see e.g. Lucy 2002, 2003). In each simulation, we follow the prop- 
agation of Monte Carlo quanta ("r-packets") representing bundles 
of X-ray photons through simply-parameterized models for accre- 
tion disk winds. During the r-packet propagation, the effects of 
Compton scattering, photoelectric absorption, bound-bound line in- 
teractions (for which the Sobolev approximation is adopted) and 
free-free absorption are simulated in detail. Interactions with matter 
(specifically bound-bound and bound-free transitions) are treated 
using the Macro Atom formalism introduced by Lucy (2002). This 
approach allows for a full treatment of line scattering, recombina- 
tion and fluorescence as required to produce realistic spectra for 
comparison with observations. The code has also been extended to 
use the method of Lucy (2002, 2003) for treating both the radia- 
tive heating and cooling of matter. In this approach, when the MC 



quanta undergo physical events in which the energy they represent 
is used to heat the plasma (e.g. by free-free absorption of photons), 
they are instantaneously converted to packets of thermal kinetic 
energy (so-called fc-packets). Assuming thermal equilibrium, these 
fc-packets are then eliminated by randomly sampling the available 
cooling processes (see Section 2.3). Depending on the cooling pro- 
cess selected, the fc-packet may be converted to an r-packet which 
is then free to propagate. The r-packets created in this way simulate 
the cooling radiation emitted by the wind. 

During each Monte Carlo simulation, the physical wind prop- 
erties (e.g. temperature, ionization conditions etc.) are held fixed. 
At the end of each such simulation, the histories of the Monte 
Carlo quanta are used to make improved estimates of the wind 
properties assuming radiative, thermal and ionization equilibrium 
and then the Monte Carlo experiment is repeated. In this way, the 
wind properties are iterated to consistency with the radiation field. 
After the iteration cycle, a final Monte Carlo simulation is per- 
formed from which viewing-angle dependent spectra are extracted 
(see Section 2.7). 



2.2 Ionization 

To provide a more complete description of the ionization state of 
the outflow, the modified nebular approximation adopted in Paper I 
has been replaced with a detailed computation of ionization and 
recombination rates that are then used to solve for the steady-state 
ionization fractions. The following sections describe the processes 
that are included. 



2.2.1 Photoionization 

Bound-free absorption involving the ejection of an outer shell elec- 
tron is included for ground states and low-energy metastable levels 
of all ions. During the Monte Carlo simulation, the packet trajecto- 
ries are used to record estimators for the photoionization rate coef- 
ficient (7) for each bound-free process in each of the computational 
grid cells of wind properties. The estimator for the rate coefficient 
of bound-free absorption from level j in a cell p is given by 
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where V p is the wind cell volume, aj is the photoabsorption cross- 
section for level j (see Section 2.9) and the summation runs over all 
Monte Carlo quanta trajectories inside cell p when the quanta have 
co-moving frame frequency v that is greater than the threshold uq. e 
is the co-moving frame energy of the quantum and ds is the trajec- 
tory length. At is the time interval represented by the Monte Carlo 
experiment. The Doppler factor correction term for transforming 
the photon path length from the observer frame to the fluid frame 
(see e.g. Lucy 2005) is neglected for computational expediency. 

At the end of a Monte Carlo simulation, the 7 estimators for 
levels of each ion are combined to make a photoionization rate es- 
timator for each ion (/) in each cell (p) defined by 



(2) 



where the summations run over all level j of the ion for which 
bound-free absorption has been included in the Monte Carlo simu- 
lation. The level populations, n^k are discussed in Section 2.3.4. 
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For L- and M-shell ions, photoabsorption involving inner shell 
electrons is also included. In the treatment of these processes, de- 
tails of the outer electron configuration are neglected and all lev- 
els of the ion are assigned the same cross section. Following inner 
shell photoabsorption, the newly made ion is generally in a highly- 
excited state that may either decay radiatively or via ejection of one 
or more Auger electrons. 

In the case of K-shell photoionization of L-shell ions, at most 
one Auger electron is expected to be produced - thus the ionization 
state may ultimately increase by either one or two. To incorporate 
this in the ionization balance, Monte Carlo estimators are used to 
obtain a rate coefficient for K-shell photoabsorption for each ion 

V>V K 

where the summation now runs over all packets having frequency 
about the K-shell edge threshold, vk and a K is the K-shell pho- 
toabsorption cross-section. To obtain separate rates for single and 
double ionization, the mean probability (p/[l]) of K-shell absorp- 
tion being followed by ejection of a single Auger electron is re- 
quired. This is computed using 

where < ... > j is a weighted mean over all states j of the ionization 
state I + 1 that are accessible by K-shell photoabsorption of ion 7, 
J2i Aji is the total autoionization rate out of level j and J2i ^iiPji 
is the sum of radiative decay rates incorporating the Sobolev es- 
cape probabilities j3ji over all downward transitions from level j. 
In practice, only states reached by K-shell photoabsorption from 
the ground configuration of ion I are included in the average and 
these are weighted with their statistical weights. The estimated rate 
at which K-shell photoabsorption in an L-shell ion increases the 
ionization state by one is then given by 

r£ fc (i- P j[i]) (5) 

and, for double ionization, 

rf>/[i] . (6) 

Data sources for Aji and are discussed in Section 2.9. 

K- and L-shell photoionization of M-shell ions, are treated in 
a similar manner, recording a Y K ^ L estimator for each such pro- 
cess in each grid cell. Owing to the complexity of the subsequent 
decays of the vacancy state (which lead to the ejection of up to 
three Auger electrons for K-shell ionization of the M-shell ions we 
consider), the Auger ejection probabilities pj[l], Pi [2] and p.r[3], 
are not computed from the atomic data set used by the code but 
taken from the yields tabulated by Kaastra & Mewe (1993). This 
approach obviates the need to follow the decay chains of the va- 
cancy states in detail but comes at the expense of assuming that 
all the relevant radiative transitions are optically thin such that the 
M-shell ion populations are unaffected by photon trapping. 



direct (Cdi) and due to collisional excitation followed by autoion- 
ization (Cea), and both radiative (a r ) and di-electronic (a d ) re- 
combination (see Section 2.9 for data sources). 



2.2.3 Ionization equilibrium 

At the end of each Monte Carlo simulation, the computed photoion- 
ization rate coefficients (V; see Section 2.2.1) are used to solve for 
an improved set of ion populations for each element, assuming ion- 
ization equilibrium and including all the processes described above. 
Since the recombination and collisional ionization terms depend on 
the adopted kinetic temperature (T e ), the ionization balance is iter- 
ated to consistency with the calculation of T e (see below) for fixed 
Monte Carlo estimators. 



2.3 Thermal balance 

In Paper I, T e was assumed to be uniform and treated as a model 
parameter. In order to address more general and realistic wind con- 
ditions, here we relax this assumption and obtain an estimate of T e 
as a function of position via a simplified treatment of the heating 
and cooling rates in the wind and the assumption of thermal equi- 
librium. 



2.3.1 Heating rates 

The treatment includes heating by Compton scattering, free-free 
absorption and bound-free absorption of the X-ray radiation field 
described by the Monte Carlo simulations. Estimators for each of 
these heating processes in each grid cell can be readily constructed 
following Lucy (2003). For Compton scattering, the heating rate in 
a grid cell (p) is given by 

where the summation runs over all packet trajectories within the 
grid cell p, a{v) is the Compton cross-section for the frequency 
(v) of the packet and f(v) is the mean energy lost per Compton 
scattering event. The free-free heating rate is 

H " =\TAt^< Kff{u)€ds ■ (8) 

Again, the summation runs over all packet trajectories in the cell. 
The free-free absorption coefficient is given by 

K ff (y) = 3£9xW s v- 3 T- 1/2 n e (l-e- h " /kT ^)J2 Z i N i cm ~ 

(9) 

where the summation runs over all ions (J) and Zi is the ion charge. 

Heating rates for photoabsorption are obtained in a simi- 
lar fashion to the ionization rate estimators in 2.2.1. For bound- 
free absorption involving outer shells, level-by-level estimators are 
recorded and used to obtain a total heating rate for each ion J, 



2.2.2 Collisional ionization and recombination processes H b / k = n j,khY k (10) 

In addition to the photoionization and Auger processes described 

above, the ionization balance includes collisional ionization, both where 
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h 5i = T^E-iM(i-?)^- <"> 



Similarly, the heating rate due to K-shell photoabsorption by ion / 
is given by 
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All other possible heating sources are neglected, including any 
contributions from outside the spectral region of the simulated ra- 
diation field or by any non-radiative processes. 



2.3.2 Cooling rates 

Cooling rates are required for both the treatment of fc-packets (see 
Section 2.6) and to obtain an estimate of the local kinetic tem- 
perature. Cooling due to radiative recombination, electron colli- 
sional excitation, free-free emission, Compton cooling by low- 
energy photons and the adiabatic expansion of the outflow are in- 
cluded in this calculation. 

The cooling rates due to spontaneous radiative recombination 
C b / k and electron collisional excitation Cf k are obtained following 
exactly Lucy (2003; his equations 31 and 33). The free-free cooling 
rate is 



C S V S = 1.426 x 10~ 27 T e 1/2 n e Z i N ' er g s cm ~ 3 s_1 ( 13 > 

i 

where we follow Lucy (2003) in setting the mean Gaunt factor to 
one. The adiabatic rate cooling rate is approximated by 



where r is the distance from the centre of the grid cell to the coordi- 
nate origin. This expression is only a crude estimate for (7 7 since it 
neglects both the geometry of the emission regions for low-energy 
photons (i.e. it assumes a centrally concentrated point source) and 
the opacity of the outflow to low-energy photons. Nevertheless, it 
provides a convenient estimate of the Compton cooling rate that in- 
troduces no additional computational demands and should be rea- 
sonably accurate, at least for the inner (hottest) regions of the wind. 



2.3.3 The kinetic temperature determination 

To obtain the kinetic temperature (T e ) for each grid cell, the radia- 
tive heating rates computed from the Monte Carlo simulation are 
compared with cooling rates to estimate the temperature at which 
thermal equilibrium, 



Hp +H* f +y^(Hj* k +H* k ) — Cp f +Cp+Cp +^2(C b j f k +Cf tk ) 
i i 

(17) 

is established. In practice, the thermal balance equation (equa- 
tion 17) is solved iteratively with the ionization balance equations 
to obtain self-consistent estimates of both the ionization fractions 
and T e for each grid cell. 



2.3.4 Excitation 

In principle, the excitation state should be obtained from the sta- 
tistical equilibrium equations. However, to do so would require 
computation and storage of estimators for all bound-bound tran- 
sitions in every cell which is prohibitive for large computational 
grids. Therefore, we adopt a very simplified treatment of excita- 
tion, namely we use the Boltzmann distribution at the local kinetic 
temperature 



Cp = 1.5n 9 fcT e V • v 



(14) 



where n 9 is the total particle density and V • v is the divergence of 
the velocity, evaluated at the midpoint of the cell. 

The cooling rate for a population of non-relativistic electrons 
due to Compton up-scattering of low-energy photons (his << kT e ) 
can be estimated from the local energy density of such photons U 1 



(15) 



where err is the Thomson cross-section (see Frank et al. (2002), 
p. 180). The cooling rate due to Compton up-scattering of photons 
is expected to be dominated by interactions with the low energy 
photons originating from the accretion disk since they should sig- 
nificantly outnumber the X-ray photons for which our detailed ra- 
diative transfer calculations are performed. An accurate treatment 
of Compton cooling would therefore require that the entire bolo- 
metric light output of the system be considered. However, since 
our primary concern is only to obtain a reasonable estimate of the 
kinetic temperature, we avoid this complication and derive an esti- 
mate for the local value of U 1 from the total bolometric luminosity 
of the source (Lboi, which is treated as an input parameter for the 
model) and an assumption of spherical geometric dilution, i.e. 
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where rii is the population of a state with statistical weight Qi and 
energy e^, relative to the grounds state (g.s.). For the wind models 
we will describe in Sections 4 and 5, our simplified treatment of 
excitation is not expected to significantly affect the Fe Ka region 
of the spectrum since the line features which form there are mostly 
associated with low excitation states of high-ionization state ma- 
terial. The soft regions ( < 1 keV) are more likely to be affected 
owing to the large numbers of transitions between excited states 
of L-shell Fe and Ni ions which contribute to the opacity at these 
photon energies. 

2.4 Monte Carlo simulations 

The radiation transport simulations are performed using the scheme 
described in Sections 3.1 to 3.3 of Paper I, modified as described 
below. 



2.4.1 Initialization of packets 

As in Paper I, an initial power-law spectrum of packet energies is 
adopted but the input spectrum now extends from 0.1 keV up to 
511 keV to allow comparison with observational constraints from 
instruments with significant effective area at relatively high photon 
energies. 
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To account for the effects of a small but non-zero angular size 
of the X-ray emission region, the packets are no longer initialized 
exactly at the coordinate origin but in a spherical region around the 
origin. The radial extent of the emission region, r er , is a parame- 
ter in the model and is generally chosen to be several gravitational 
radii, as appropriate for an X-ray emission region of size compara- 
ble to the innermost radii of an accretion disk. 

2.4.2 Propagation of packets 

Compton scattering, bound-free absorption and line absorption are 
treated exactly as described in Paper I, except that bound-free ab- 
sorption by a few low-lying metastable levels is included for most 
ions. 

Free-free opacity from all ions is included using equation 9. 
Following free-free absorption, packets are always converted to 
thermal energy (fc-packets; see Lucy 2003). 

K-shell photoionization is included as an opacity source for all 
L-shell and M-shell ions. Details of the outer electron configuration 
are neglected such that all states of an ion contribute to the same 
opacity term. When K-shell photoabsorption occurs in the simula- 
tions, the packet either activates a macro atom or is converted to a 
fc-packet - this is exactly analogous to the treatment of outer-shell 
bound-free processes which are discussed in detail by Lucy (2003). 
Here, the probability of fc-packet conversion is simply given by 

Pk = 1 - vk/v (19) 

where v is the absorption frequency and vk is the relevant K-shell 
absorption edge frequency. If the outcome of K-shell photoioniza- 
tion is activation of a macro atom, the state activated must have a 
K-shell vacancy in its configuration and will generally be highly ex- 
cited 1 . Which of the K-shell vacancy states in the newly formed ion 
is activated is chosen randomly by sampling the statistical weights 
of the viable states in the atomic data set. 

2.4.3 Macro Atom processes 

The macro atom treatment used in Paper I has been extended to in- 
corporate autoionization, bound-free transitions from excited states 
and K-shell photoionization. 

Autoionization allows activated macro atom states to de- 
excite by conversion to a fc-packet and to make internal transition 
to higher ionization states. Following the argumentation of Lucy 
(2002), the deactivation probability (macro atom — > fc-packet) 

p D = NiA a (ei - e f ) (20) 

where A a is the rate coefficient ([s -1 ]) for autoionization from ini- 
tial state i to final state /, e» and e j are the energies (excitation plus 
ionization) of the states, and Ni is the normalization factor. Simi- 
larly, the probability of making an internal macro atom jump i — ¥ f 
is 

P L f =NiA a ef . (21) 
This formulation requires that the autoionization process connects 

1 Conversely, following activation by outer-shell bound-free absorption, 
macro atoms are always assumed to be in the ground state of the newly 
formed ion. 



two well-defined states in the atomic data set. In our calculations, 
this is done for L-shell ions using level-to-level autoionization rates 
(see Section 2.9 for data sources). However, for the M-shell ions we 
do not consider the target level in detail but assume that the energy 
flow in the autoionization process is dominated by the energy car- 
ried by the Auger electron - i.e. we assume that e/ < < e» such that 
p J /p D « 0. This has the advantage that we do not require atomic 
level-to-level autoionization rates but only total Auger-widths for 
vacancy states, as are available from the literature (see Section 2.9). 

For L-shell ions di-electronic recombination, is also included 
in the macro atom scheme. Di-electronic rate coefficients, ad, are 
obtained from A a in the usual manner and used to formulate an 
internal transition probability between the states 

p/_+/ = N f a d e f n e (22) 

where n e is the local electron number density. Note that, as for e.g. 
bound-free absorption, there is no macro-atom deactivation proba- 
bility associated with qj. 

The 7 estimators described in Section 2.2. 1 are used to com- 
pute bound-free macro-atom internal transition probabilities fol- 
lowing Lucy (2003). In all such cases, it is assumed that the rel- 
evant upper state is the ground state of the target ion. Similarly, the 
■y K estimators are used to include internal transitions associated 
with inner shell ionization. These are applied to all energy levels 
of the absorbing ion lying below the first ionization potential and 
connect to K- vacancy states of the next higher ionization state. The 
choice of which K-vacancy state to activate following an internal 
transition is made by randomly sampling their statistical weights. 
L-shell ionization of M-shell ions is not included in the macro atom 
treatment. 

2.5 Accretion disk 

We assume that the black hole is surrounded by an optically thick 
accretion disk that lies in the xy-plane and extends from an inner ra- 
dius Td outwards throughout the simulation domain. r d is treated as 
a model parameter but set equal to 6 r 9 (where r 3 = GM bh /c 2 ) for 
all the simulations performed here. During the Monte Carlo simula- 
tions, we also assume that all photon packets that strike the disk are 
lost to the X-ray regime and they are removed from the simulation. 

For our particular choice of source geometry, the disk sub- 
tends a modest solid angle as seen by the X-ray source (~ 2.8 sr), 
meaning that a fraction of the primary X-ray photons strike the in- 
ner regions of the disk directly, without interaction in the wind. 2 
The total flux striking the disk, however, is greater than this since 
the wind effectively causes irradiation of the outer parts of the disk 
via scattering and emission in the wind. Physically, these photons 
will interact with the disk and may lead to an additional compo- 
nent of disk reflection in the X-ray spectrum. We neglect any such 
component here but it will be investigated in a future study. 

2.6 Thermal emission and fc-packets 

In paper I, it was assumed that energy packets converted to fc- 
packets would be lost to the hard X-ray band. However, since cool- 
ing rates are now computed in order to impose thermal equilibrium 

2 This fraction, however, is not physically meaningful in our model since 
it is merely determined by our simple choice for the geometry of the X-ray 
emitting region. 
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(Section 2.3.2), these are now used to re-emit the packet, thus sim- 
ulating the thermal emission of the outflow. 

When conversion to a fc-packet occurs, the cooling terms de- 
scribed in Section 2.3.2 are first randomly sampled to choose a 
cooling process. If C a is chosen, the energy is assumed to be lost to 
the radiation field. The collisional cooling processes (C cl ) lead to 
activation of a macro atom state, the particular state being chosen 
by sampling the terms contributing to C cl (see Lucy 2003). The full 
macro atom machinery then determines the outcome of this event. 

Choosing C ib or leads to conversion to photon packets. 
Both processes are assumed to emit isotropically in the co-moving 
frame and the emission frequency is selected using equations 26 
and 41 of Lucy (2003). Note that this emission rule for bound- 
free processes is not exact here since it assumes a simplified cross- 
section shape but it is significantly less computationally demand- 
ing that sampling the real emissivity and is adequate for current 
purposes. 

2.7 Extraction of spectra 

Although spectra may be obtained by directly binning the emergent 
packets by frequency and direction of travel, as done in Paper I, 
this has the disadvantage that very large numbers of quanta must 
be simulated to suppress Monte Carlo noise, particularly when the 
spectrum is expected to depend strongly on the observer line-of- 
sight. As discussed by e.g. Lucy (1999), that approach is far from 
optimal and does not make full use of the packet trajectories during 
the Monte Carlo simulation. Ideally, volume-based Monte Carlo es- 
timators should be used to obtain emissivities throughout the vol- 
ume of the simulation and used to obtain spectra via a formal so- 
lution of the radiative transfer equation (see discussion e.g. Lucy 
1999). However, this approach is very demanding of memory re- 
sources: for our atomic data set (with 10 5 atomic processes) and 
computational grids (lOOx 100 wind cells), it ideally requires the 
storage of ~ 10 9 floating-point estimators which is prohibitive. 
Therefore we extract spectra using a approach similar to that im- 
plemented in the disk-wind simulations of Long & Knigge (2002) 
which, although less efficient, does not make significant demands 
on memory consumption. 

Viewing-angle dependent spectra are only extracted during the 
final Monte Carlo simulation, once the ionization and thermal prop- 
erties of the outflow have been iterated to consistency with the radi- 
ation field. The set of lines-of-sight for which spectra are required 
is pre-specified and then every physical event (meaning emission 
or scattering of a photon-packet) that occurs in the Monte Carlo 
simulations is used to compute a contribution to each of the spectra 
in the following manner. 

When a physical event takes place in the Monte Carlo simula- 
tion, the propagation of the packet is temporarily suspended and the 
probability per unit solid angle of the the associated photon-packet 
having been scattered into each of the requested lines-of-sight is 
computed, dP(n)/dfi where the vector n identifies the line-of- 
sight. Depending on the physical process responsible for the packet 
emission/scattering, dP(n)/dQ may depend on previous proper- 
ties of the packet (e.g. its incoming direction in the case of Comp- 
ton scattering) or on the local outflow properties (e.g. the Sobolev 
escape probability in the case of line emission). We then integrate 
the total optical depth, r (n), from the position of the physical event 
to the edge of the wind in the direction n. r(n) is computed ex- 
actly as in the Monte Carlo simulation, incorporating contributions 
from all physical processes included in the code and allowing for 
Doppler shifts to bring the photons into resonance with spectral 



Table 1. Elements and ions that are included in the radiative transfer calcu- 
lations. 



Element 


Ions 


Element 


Ions 


C 


IV - VII 


S 


VIII - XVII 


N 


V - VIII 


Ar 


X - XIX 


O 


VI - IX 


Ca 


XII - XXI 


Ne 


VIII - XI 


Fe 


X - XXVII 


Mg 


IV - XIII 


Ni 


XII - XXIX 


Si 


VI - XV 







lines. With dP(n)/dSl and r(n) computed, we then add an en- 
ergy contribution from the event to the spectrum associated with 
direction n of 



where e r f is the rest-frame energy of the packet and the contribu- 
tion is added to the frequency bin associated with the rest-frame 
frequency of the packet. This calculation is repeated for each of 
the lines-of-sight for which spectra have been requested and then 
the Monte Carlo simulation proceeds until the next physical event 
occurs. In this way, we obtain high signal-to-noise spectra without 
resorting to coarse angular binning of the emergent spectra or re- 
quiring the storage of prohibitively large numbers of Monte Carlo 
estimators. 

2.8 Rescaling of model parameters 

Although our numerical simulations are performed adopting a par- 
ticular value for the mass of the central black hole (Mbh), the spec- 
tra obtained are applicable for other Mbh -values following a rescal- 
ing of the dimensional wind parameters. Specifically, the ionization 
state, optical depths and velocity-law in the wind are preserved un- 
der a global rescaling where all the model luminosity parameters 
(Lboi, Lx) scale with the Eddington luminosity (Z/Edd), all lengths 
(r m in, r max , d, R v , Td and r em ) are scaled to the gravitational radius 
(r 9 ) and the mass-loss rate (M) is scaled to the Eddington accretion 
rate (MEdd). 

2.9 Atomic data 

Table 1 lists the elements and ionization stages for which atomic 
models are included in the radiative transfer calculations. Com- 
pared to Paper I, the L-shell ions of the astrophysically abundant 
intermediate-mass and iron group elements have been added (Mg, 
Si, S, Ar, Ca, Fe and Ni). The highest M-shell ions (down to the 
Cl-like ion) have also been added for Fe and Ni. The solar element 
abundances of Asplund et al. (2005) are adopted. 

In all cases, inner-shell photoionization cross-sections were 
taken from the fits by Verner & Yakovlev (1995). We note that, al- 
though convenient, those data do not account for smoothing of the 
K edges by broad resonances below threshold (see Palmeri et al. 
2002; Kallman et al. 2004). Depending on the ionization state, this 
effect could be critical for quantitative analysis of the region around 
the Fe K edge but it is not important for the He- and H-like ions 
which dominant in most of the wind models we describe in Sec- 
tions 4 and 5. For outer-shell photoionization from ground config- 
urations (and also low-lying metastable states), cross-sections were 
taken from Verner et al. (1996) and a simple hydrogenic approxi- 
mation was used for cross-sections from more highly excited states. 
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Atomic energy levels, bound-bound A-values and electron 
collision strengths were extracted from the CHIANTI atomic 
database, version 5 (Dere et al. 1997; Landi et al. 2006) for all 
the included L- and M-shell ions of C to Fe. For the M-shell ions 
of Fe, the data were limited to those of the configurations included 
in the study by Mendoza et al. (2004). In all cases, the T e depen- 
dence of the electron collision strength was neglected and the low- 
temperature limit adopted. 

Autoionization rates for K- vacancy states in L-shell ions of Fe 
were taken from Palmeri et al. (2003) and Auger widths for the K- 
vacancy states of the M-shell Fe ions from Mendoza et al. (2004). 
The total autoionization rates for the K- vacancy states of Mg, Si, S, 
Ar and Ca were taken from Palmeri et al. (2008) and the branching 
ratios of the autoionization process into the specific final states of 
the target ion were assumed to be the same for the iso-electronic Fe 
ion. 

For the included L- and M-shell ions of Ni (XII - XIX), energy 
levels, Einstein A-values and total Auger widths were taken from 
Palmeri et al. (2008). Again, the L-shell autoionization branching 
ratios were assumed to be the same as in the iso-electronic Fe ion. 

Radiative recombination rates are taken from Gu (2003b), 
where available, or else from Verner & Ferland (1996) (H-, He- 
and Li-like ions of C, N, O and Ne and Na-like ions of Fe and Ni), 
Arnaud & Raymond (1992) (Fe XVI) or Shull & van Steenberg 
(1982b,a), otherwise. Di-electronic rates were obtained from Gu 
(2003a), supplemented by those from Arnaud & Raymond (1992) 
(M-shell Fe ions), Shull & van Steenberg (1982b) (M-shell Ni ions) 
and Arnaud & Rothenftug (1985) (He- and Li-like C, N, O and 
Ne). Collisional ionization rate coefficients (Cdi and Cea) are 
obtained from Arnaud & Raymond (1992) for Fe and Arnaud & 
Rothenflug (1985) for all other elements. 

2.10 Wind models 

The modification described above mean that the code is now ca- 
pable of simulating all the necessary atomic physics to compute 
X-ray spectra for wind models and several such calculations will 
be presented in the Sections below. 

For the moment, we continue to work with the simply- 
parameterized, smooth outflow models exactly as described in Pa- 
per I. This approach allows us to investigate outflow signatures and 
their sensitivity to the major wind parameters in a relatively simple 
manner. In the future, however, we plan to extend our studies to 
models that go beyond the smooth, steady-state flow prescriptions 
we adopt here. Numerical simulations clearly suggest that outflows 
are likely to have complex structure that will affect the observed 
spectra (see e.g. Proga & Kallman 2004; Schurch et al. 2009) and 
time variation in the outflow properties may well have a role in ex- 
plaining the observed time dependence of some spectral features. 



3 TEST OF IONIZATION BALANCE 

The most significant improvements in the physics of the current 
code compared to that described in Paper I are the extension to 
lower ionization states and the substantially more sophisticated 
treatment of the ionization balance. To test these improvements, we 
have made a simple comparison with the ionization balance in the 
well-known photoionization radiative transfer code CLOUDY (ver- 
sion 07.02, last described by Ferland et al. 1998). For the test, we 
used a wind model with very low column density such that it would 
be optically thin and every cell would be illuminated with the pure 
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Figure 1. Ionization fractions as a function of ionization parameter for Fe 
XVI - XXVII computed for optically thin conditions with CLOUDY (solid 
lines) and our Monte Carlo code (diamond symbols). In both calculations, a 
power-law continuum with photon index T = 1.2 and a kinetic temperature 
of T e = 10 6 K were adopted. 

power-law continuum from the X-ray point source. A primary pho- 
ton power-law index of V = 1.2 was used for the test and the ki- 
netic temperature was forced to be 10 6 K everywhere. Ionization 
fractions were computed using the full machinery of our code and 
the results for Fe ions are plotted in Figure 1 versus an ionization 
parameter defined as 

C=^^ (24) 

where Lo.i-so is the source luminosity between 0.1 and 50 keV. 
The Fe ionization fractions obtained with CLOUDY for a low 
column-density spherical shell are over-plotted for comparison. 

The agreement between the ionization calculations in the two 
codes is very good, better than 10 per cent for most of the K- and 
L-shell ions of iron. This precision is comparable to the accuracy 
of the source atomic data. Our calculations become less reliable for 
ionization parameters much below f ~ 100 ergs cm s _1 since our 
current atomic data set does not contain all the necessary ions (in 
our test calculation, Fe IX - which we neglect - has a significant 
role in the CLOUDY ionization balance below £ ~ 50 ergs cm s _1 ). 



4 EXAMPLE CALCULATIONS 
4.1 Model parameters 

Here we present details of two outflow models computed with the 
improved code. The first, Model A has the same outflow parameters 
as the example model described in Paper I. However, the physical 
properties of the flow differ slightly from Paper I since the kinetic 
temperature is no longer a parameter but is determined as described 
in Section 2.3.3. Also, since the new code version allows for a finite 
spatial extent of the primary X-ray emission region and occultation 
of the scattered/reprocessed light (see Section 2.5), the relative am- 
plitudes of the direct and scattered/reprocessed components of the 
spectrum are altered (see Section 4.3). 

Our second model, Model B, has parameters chosen to illus- 
trate the effects of lower ionization state material in the wind as 
can now be treated with the improved code version. We adopted 
a black-hole mass of 10 7 Mq and an outflow that is wider than 
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those considered in Paper I, extending across a factor of three range 
of launching radius on the disk (50 - 150 r g ). Since the ioniz- 
ing source is centrally concentrated this means that the ionization 
gradient across the wind (which was already present in Paper I) 
now extends to even lower ionization states around the outer edge. 
We have also chosen a more gradually accelerated outflow with 
lower terminal velocity. Specifically, we adopt a larger value of 
the velocity-law acceleration length parameter R v (which speci- 
fies how far downstream the outflow reaches one half of the ter- 
minal speed; see equation 1 of Paper I) and a smaller value of /„ 
(which relates the terminal speed to the escape speed at the base of 
a streamline; see equation 2 of Paper I). The X-ray source param- 
eters and the wind mass-loss rate for Model B were motivated by 
previous studies of the bright quasar PG121 1+143 with which more 
detailed comparisons are made in Section 5. The primary X-ray 
source luminosity and the total bolometric luminosity were chosen 
to give ratios to the Eddington luminosity (Lndd) similar to those 
inferred for PG121 1+143: adopting a typical 2-10 keV luminos- 
ity of 10 44 ergs s _1 and bolometric luminosity of 4 x 10 45 ergs s _1 
(Pounds et al. 2003) yields L x /Lba& ~ 0.02 and L m /Lem ~ 0.8 
for a black hole mass ~ 4 x 10 7 M Q (Kaspi et al. 2000) and ac- 
cretion efficiency of ~ 0.06. The power-law index (T ~ 1.8) and 
mass-loss rate (M ~ Mem) were also chosen to be roughly appro- 
priate for PG121 1+143, as motivated by the discussion of Pounds 
et al. (2003). 

The complete set of model parameters is given in Table 2; def- 
initions of the geometry and velocity-law parameters are given in 
Paper I. For both models, we placed the innermost edge for the ac- 
cretion disk at the last stable orbit for a Schwarzschild black hole 
(6r 9 ) and assumed a comparable size for the region of X-ray emis- 
sion (again 6r 9 ). 

4.2 Thermal and ionization structure 

Figure 2 shows the distribution of kinetic temperature and Fe ion- 
ization states in the example models. For Model A, the spatial vari- 
ation of the ionization state is qualitatively similar to that discussed 
in Paper I: the inner edge of the wind is almost fully ionized and 
gradients of decreasing ionization occur both along the outflow and 
across the base of the flow. Model B is similar but, owing to the 
larger radial extent of the flow launching region in this model, the 
gradient across the flow is particularly well-developed such that a 
substantial region on the outside of the wind is dominated by L- 
shell Fe ions. As expected, the regions of the wind that are most 
strongly ionized are also most strongly heated by the X-ray source 
leading to significant kinetic temperature variations across the flow 
(right panels of Figure 2). 

4.3 Computed spectra 

For both the example models, we computed emergent spectra for 
twenty lines-of-sight uniformly sampling the orientation (0 < fi < 
1, where 6 = cos -1 /i is the angle between the line-of-sight and 
the rotation axis). Figures 3 and 4 show the computed spectra for 
five of these lines-of-sight for Models A and B, respectively. The 
orientations of the chosen lines-of-sight are indicated by the dashed 
lines in Figure 2. As in Paper I, the propagation histories of the 
Monte Carlo quanta are used to divide the emergent spectrum into 
a direct component (representing photons that reached the observer 
without any interactions) and a scattered/reprocessed component 
derived from the quanta that underwent at least one interaction in 



the wind. These contribution are separately plotted in the left panels 
of Figures 3 and 4. The right panels show the 2 - 10 keV region of 
the spectra in greater detail. 

The first line of sight shown for both models (top panels of 
Figures 3 and 4) corresponds to an observer at sufficiently low in- 
clination that no portion of the wind obscures the primary X-ray 
source. Thus the direct component of radiation is unaffected by the 
wind (red line in top left panels). However, since the outflow scat- 
ters light from other lines-of-sight and produces its own thermal 
emission, the direct light is supplemented by a significant compo- 
nent of scattered/reprocessed radiation. As one would expect, this 
scattered/reprocessed spectrum is qualitatively similar to that ob- 
tained from standard disk reflection models (e.g. Ross & Fabian 
2005). In particular, a moderately strong Fe Ka emission line ap- 
pears in both models. This is predominantly formed by Fe XXV and 
XXVI since most of the wind reflection seen from this line-of-sight 
occurs in the highly-ionized inner surface of the outflow. The emis- 
sion line profile is Doppler broadened and, owing to the substantial 
outflow velocities, develops an electron scattering wing (Auer & 
van Blerkom 1972; Laurent & Titarchuk 2007) that skews the pro- 
file to the red (for both models, it extends down to ~ 5 keV for this 
line-of-sight). At softer energies, the spectrum of scattered Monte 
Carlo quanta also contains a variety of other emission features. 
These are considerably stronger in Model B owing to the greater 
range of ionization conditions in the wind. In this model, weak but 
distinct Ka emission lines of O, Si and S (~ 0.65, 2.0 and 2.6 keV, 
respectively) are present together with forests of blended lines from 
the L-shell ions of Fe. The emission profiles of these features are 
broadened by the same mechanisms which affect the Fe K line lead- 
ing to rather complex spectra. The scattered/reprocessed spectrum 
also adds a distinct "Compton hump" to the total spectrum, caus- 
ing it to have a broad peak around 20 - 30 keV; at even harder 
energies the scattered/reprocessed spectrum bends down (owing to 
the energy dependence of the Compton cross-section) such that the 
flat, primary X-ray spectrum becomes increasingly dominant above 
about 50 keV. For this inclination angle, the spectrum computed for 
Model A is very similar to that obtained for the same outflow pa- 
rameters using the code version in Paper I (see upper right panel of 
Figure 3). The main difference is a modest offset in normalization 
that arises from the reduced amplitude of the scattered component 
of radiation owing to occultation by the accretion disk which is in- 
cluded in the new calculations (see Section 2.5). 

The second line of sight considered (p — 0.675) passes 
through the upper-most layers of the outflow models which are 
hottest and most highly-ionized. The remaining three lines-of-sight 
(jj, = 0.575, 0.475 and 0.275) pass through increasingly denser and 
cooler parts of the wind such that the spectra become more com- 
plex. As in Paper I, lines of sight through the wind show weak, 
narrow, blueshifted absorption features associated with highly ion- 
ized material, most importantly the Ka lines of Fe XXV and XXVI. 
In Model A, the spectra for /i — 0.675 (and lower) are charac- 
teristically similar to that of the example model in Paper I (see 
comparison in Figure 3) but there are some differences in detail. 
In particular, with the improved code version the amplitude of the 
component of scattered/reprocessed light is reduced by disk occul- 
tation, the mean ionization state is slightly lower and the absorption 
features are sharper. This last effect arises because of the angular 
binning of MC quanta used to extract the spectra in Paper I: since 
the line-shifts depend on inclination, spectra obtained via angular 
binning will tend to have slightly smeared line profiles. 

In Model B, the = 0.675 line of sight has an integrated hy- 
drogen column density of Nh ~ 2 x 10 23 cm~ 2 , only a factor of 
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Table 2. Inputs parameters for the example model. 



Parameter 



Model A 



Model B 



mass of central object, Mj,h 

bolometric source luminosity, Lb„i 

source luminosity (2-10 keV), Lx 

source power-law photon index, V 

range of source photon energies in simulation 

size of emission region, r er 

inner radius of disk,r d 

inner launch radius, 

outer launch radius, r m ax 

distance to wind focus, d 

terminal velocity parameter, f v 

velocity scale length, R v 

velocity exponent, j3 

launch velocity, vo 

wind mass-loss rate, M 

mass-loss exponent, k 

outer radius of simulation grid 

3D Cartesian RT grid cells 

2D wind grid zones 



4.3 X 10 b M 

2.5 X 10 44 ergs s" 1 (~ 0.6L Edd ) 
10 43 ergs s- 1 (~ 0.02L Edd ) 
2.38 

0.1-511 keV 

6r 9 = 3.8 X 10 12 cm 
6r = 3.8 X 10 12 cm 



6.4 X 10 13 cm 



100 r 9 = 



' min 

1.0 

An in 

1.0 

o.o 

0.1 M yr" 1 (~ 0.6M Edd ) a 
-1.0 

5 X 10 16 cm 
180 x 180 x 180 
100 x 100 



10 7 Mq 

10 45 ergss" 1 ^ 0.8L Edd ) 

2.5 X 10 43 ergs s" 1 (~ 0.02L Edd ) 

1.8 

0.1-511 keV 



6r 9 : 
6r 9 : 
50 r 9 

3r, 



■ 8.9 x 10 12 cm 
: 8.9 x 10 12 cm 
= 7.4 x 10 13 cm 



111111 

min 



0.5 

3^m in 

1.0 
0.0 

0.38 M yr" 1 (~ M Edd )° 
-1.0 

5 X 10 16 cm 
180 x 180 x 180 
100 x 100 



The Eddington accretion rate, Af Edd was computed assuming a radiative efficiency of six percent. 
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Figure 2. Distribution of mean Fe ionization state (left) and kinetic temperature (right) for Models A (top) and B (bottom). The dashed lines indicate the five 
lines-of-sight for which spectra are shown in Figures 3 and 4. 
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Figure 3. Spectra for five observer lines-of-sight computed from Model A (from top to bottom, = 0.875, 0.675, 0.575, 0.475 and 0.275 where 8 = cos i± 
is the angle of the line-of-sight relative to the polar axis). The left panels show the spectrum from 0.1 - 200 keV while the right panels show the 2 - 10 keV 
region in detail for the same lines-of-sight. Note that the abscissa is plotted logarithmically in the left panels but linearly in the right panels. In the left panels, 
the total spectrum is shown in black, the spectrum of direct photons in red and the scattered/reprocessed spectrum in blue. In the right panels, the total spectrum 
computed with the new code version is plotted in black while the spectra from the example model in Paper I are shown in green, for comparison. The vertical 
lines indicate the mean rest energies of the Fe XXV/XXVI Kct transitions (~6.7/6.97 keV). All the spectra are normalised to the incident X-ray spectrum which 
has a photon power-law index, T=2.38. Monte Carlo noise from the simulations is present in all the spectra and is responsible for the small-scale fluctuations. 
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Figure 4. As Figure 3 but showing spectra from Model B. No comparison spectra are shown in the right panels. All the spectra are normalised to the incident 
X-ray spectrum which has a photon power-law index, r=1.8. 



a few smaller than that suggested for the most highly ionized ma- 
terial in the spectrum of PG121 1+143 (~ 5 x 10 23 cm -2 , Pounds 
et al. 2003). However, the material along this line of sight is hot and 
almost fully ionized (see Figure 2) such that Compton scattering 
is the only important opacity source. This reduces the flux of the 
direct component of the radiation and causes it to curve upwards 
at high energies, a consequence of the photon-energy dependence 
of the Compton cross section. Although the scattered/reprocessed 
component of the light still shows a significant Compton hump, 



this upwards curvature of the direct component partially cancels it 
out leaving only a modest bump in the total spectrum. The only 
significant absorption feature imprinted by the wind for this line 
of sight is from the trace population of Fe XXVI and only the Ka 
line is strong enough to be plausibly detectable with current instru- 
ments. The emission features in the scattered/reprocessed spectrum 
are qualitatively very similar to those for the lower inclination angle 
described above except that the features are slightly broader thanks 
to the greater line-of-sight components of both the outflow and ro- 
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tational velocity fields for this orientation. In addition, the Fe Ka 
line shows two distinct peaks, one close to the rest Ka energy of 
Fe XXVI and the other to that of Fe XXV and the high L-shell ions. 

The third to fifth lines-of-sight for which spectra are shown 
for Model B (Figure 4) have sufficiently high column densities 
(N H ~ 10 24 , ~ 4 x 10 24 and ~ 10 25 cm~ 2 , respectively) that 
the direct component of radiation is almost completely blocked ex- 
cept at rather high energies ( > 30 keV). For the highest inclination 
angle considered (p = 0.275), the reprocessed component of light 
dominates even well above 100 keV such that the total spectrum of 
the wind shows a downturn at these energies, a property that can be 
constrained by high-energy observations (see discussion by Fabian 
et al. 1995). However, we stress that for the lower inclination angles 
in this model, and indeed all the angles considered for Model A, the 
direct component rises sufficiently rapidly with energy that there is 
still significant flux emerging above 200 keV. 

For the fi = 0.575, 0.475 and 0.275 lines-of-sight, the 
Model B spectrum is progressively more strongly affected by ab- 
sorption in the wind. The Fe Ka absorption line becomes deeper 
and gradually switches from being dominated by Fe XXVI to 
Fe XXV at higher inclination angles. The associated Ka emission 
line remains strong and its apparent red wing becomes increasingly 
extended for high inclination angles. As expected, the Fe Ka ab- 
sorption is also accompanied by blueshifted absorption by Fe K/3 
(~ 9.0 keV) and Ni Ka (Ni XXVII at ~ 8.5 keV and XXVIII 
at ~ 8.8 keV). There is also significant absorption around the 
blueshifted Fe K edge which, although somewhat smeared out by 
the Doppler shifts, causes the distinct drop in flux at around 10 keV. 

At softer energies, blueshifted absorption features associated 
with the abundant light elements become increasingly prominent as 
the inclination angle is increased. For an inclination of pi = 0.475 
(fourth row in Figure 2), weak narrow Ka absorption lines of S XVI 
(~ 2.8 keV), Si XIV (~ 2.2 keV), Mg XII (~ 1.6 keV), Ne X 
(~ 1.1 keV) and O VIII (~ 0.71 keV) are all present but their 
equivalent widths are small, generally less than a few eV. Weak 
Fe XXIV lines also appear at ~ 1.2 and 1.3 keV due to blueshifted 
transitions between its ground 2p° configuration and states of the 
excited 2p 4 ( 3 P)3s configuration. At the higher inclination angle of 
fj, — 0.275 (fifth row in Figure 4), the absorption spectrum is even 
more complex. The Ka transitions mentioned above are all still 
present and generally somewhat stronger. Very weak Ka lines due 
to Ar and Ca also manifest (~ 3.6 and 4.4 keV, respectively), as do 
the K/3 lines of the H-like ions of Mg, Si and S. Although O VIII 
K/3 absorption is also significant in the simulation, this feature is 
blended with comparably strong lines of Fe XVIII and XIX. Deep 
absorption occurs between 1.0 and 1.3 keV owing to the forest of 
lines arising in transitions from the low-lying n — 2 configurations 
of Fe and Ni L-shell ions. Furthermore, there is also significant 
bound-free absorption by the ground states of the K-shell ions of 
O, Mg, Si and S and of the L-shell ions of Fe that causes a drop in 
the flux above ~ 1.5 keV. 



5 PRELIMINARY COMPARISON WITH PG1211+143 

Thanks to the extension of the capabilities of the code to deal with 
lower ionization conditions, we can now use the models to compare 
much more realistic theoretical spectra with observations to estab- 
lishing how readily outflow models can explain observed spectro- 
scopic features and, if real, what the physical properties of these 
flows might be. For a fully detailed comparison a large grid of the- 
oretical models covering the physically interesting parameter space 



must be computed. This will be an important next step in our stud- 
ies but we begin here by making a preliminary comparison with the 
observed spectrum of the bright quasar PG121 1+143. 

5.1 Observational data 

Amongst the first direct evidence for highly ionized massive out- 
flows from AGN was a report by Pounds et al. (2003, 2005) based 
on the analysis of XMM-Newton X-ray spectra of PG121 1 + 143. 
They reported the detection of several strong absorption features 
which they associated with blueshifted Ka transitions of C, N, O, 
Ne, Mg, S and Fe. To quantify the absorption features in the (1 
-10 keV) European Photon Imaging Camera (EPIC) pn spectrum 
they fit an XSTAR (Kallman & Bautista 2001) absorption model 
and found that a high column density (Nh ~ 5 x 10 23 cm -2 ) and 
a high ionization parameter (log(£/ergs cm s _1 ) ~ 3.4) were re- 
quired to explain both the claimed Fe and S Ka absorption lines as 
originating in a rapid (v ~ 0.08c) outflow. Based on the proposed 
identification of the lower ionization Mg line in the EPIC pn spec- 
trum and lines attributed to the lighter elements in their contempo- 
raneous Reflection Grating Spectrometer (RGS) spectra, they also 
postulated that additional absorbing components must be present in 
the outflow with significantly lower column density and ionization 
parameter (although comparable outflow speeds). Broad emission 
was also identified in the EPIC pn spectrum, extending from ~ 3 
to 7 keV. Although this emission feature can be fit as an extreme 
relativistic Fe Ka emission line originating from reflection by the 
inner most regions of an ionized accretion disk (Pounds et al. 2003), 
it has more recently been suggested that it may be P Cygni emis- 
sion physically associated with the blue-shifted Ka absorption (e.g. 
Pounds & Reeves 2009) as would be expected from a wind with a 
large covering fraction. 

Subsequent re-observations of PG121 1+143 with XMM- 
Newton were made in 2004 (Pounds & Reeves 2007) and 2007 
(Pounds & Reeves 2009). These showed that, although the spec- 
tral properties of PG121 1+143 are significantly time- variable, the 
fast outflow of highly ionized gas is persistent and confirmed that it 
is likely to be an energetically important component of the system. 

Given the complexity of the available data for PG121 1+143 
and the wide range of spectroscopic features identified in the 2001 
XMM data first reported by Pounds et al. (2003), this object pro- 
vides an excellent point of comparison with our synthetic spec- 
tra for outflows. In particular, while previous efforts to quantify 
the spectra have generally adopted distinct fit components for the 
claimed absorption lines and the apparently broad emission red- 
ward of ~7 keV, our wind models allow us to explore the relation- 
ship between these features with a set of physical models. 

The X-ray spectrum of PG121 1+143 shows a strong excess 
for energies < 1 keV relative to an extrapolation of a power-law 
fit to the harder energy spectrum (see e.g. Pounds et al. 2003). The 
origin of soft excesses in AGN spectra remains unclear: although it 
has been suggested that they may be a consequence of absorption 
by relatively low ionization material and data showing strong soft 
excesses (including those of PG121 1+143; see Pounds & Reeves 
2007, 2009) have been successfully fit assuming an absorption- 
dominated origin, the geometry and origin of the absorbing ma- 
terial is not well-constrained. In view of this uncertainty, we re- 
strict our fitting to the 2-10 keV spectrum of PG121 1+143 to 
avoid biasing our fit by a choice of particular model for the ori- 
gin of the soft excess. However, we will qualitatively discuss the 
lower-energy absorption features predicted by the model compared 
to those reported in the PG121 1+143 spectra. 
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In the following, we will make quantitative comparisons with 
the data set compiled by Pounds & Reeves (2009) by stacking XMM 
EPIC pn data from 2001, 2004 and 2007. We have rebinned the data 
to the energy-dependent half-width at half-maximum (HWHM) of 
the EPIC pn instrument. 

5.2 Choice of model parameters 

We began by making a comparison of the stacked 2-10 keV data 
of PG121 1+143 with our example spectra (described in Section 4). 
As discussed above, Model B already indicated that many of the 
discreet features identified in the spectrum of PG121 1+143 can 
qualitatively be reproduced in our theoretical spectra. To produce a 
strong, significantly blueshifted Fe XXV/XXVI absorption feature, 
the example models suggest that a line-of-sight passing through 
some portion of the inner, highly ionized edge of the flow is re- 
quired (fx ~ 0.4 - 0.6 for the example model geometric parame- 
ters). Although Model B already produces strong Ka absorption, 
the blueshift is slightly too small compared to that suggested for 
PG121 1+143 (Pounds et al. 2003, 2005; Pounds & Page 2006). 
Therefore, we first increased the outflow velocity parameter (/„) 
from 1/2 to 2/3, a value intermediate between those used for the 
two example models. 

Upon quantitative comparison, the Model B spectra also dif- 
fer in detail from the observations of PG 121 1+143 in that the lower 
ionization features are somewhat too weak for the appropriate in- 
clination angles (e.g. the equivalent width of the S XVI Lya line is 
only < 5eV in the model while Pounds et al. (2005) and Pounds 
& Page (2006) report equivalent widths of ~ 40 and ~ 24eV, re- 
spectively). Thus, models with lower typical ionization state are 
desirable. In principle, a full grid of models exploring all the out- 
flow parameters should be explored to find the best match to the 
observed spectrum. Here, however, we restrict ourselves to varying 
parameters that are expected to directly affect the mean ionization 
state. We therefore created a small grid of six models that explores 
the effects of higher wind mass-loss rates (M) and longer velocity- 
law scale lengths (R v ) while keeping all other parameters the same 
as for Model B. Relatively modest numbers of Monte Carlo quanta 
were used for these simulations in order to reduce the computation 
demands (only ~ 10 6 packets were used for the calculation of the 
spectra). 

The quantitative effect of M on the models was discussed in 
Paper I and that discussion remains applicable here. In particular, 
increasing M leads to higher densities and therefore lower ioniza- 
tion states for a fixed X-ray source luminosity. In general, it also 
leads to stronger absorption features thanks to the increased col- 
umn density for a fixed line of sight. Our grid of models explored 
M-values of 0.38, 0.76 and 1.52 M Q yr _1 corresponding to ~ 1, 2 
and 4 Me<m for our adopted value of Mbh ~ 10 7 Mq. 

The second parameter varied, R v , affects the spectrum for two 
distinct reasons. First, a larger value of R v means that the flow 
accelerates more gradually such that the material in the accelerating 
region is more dense for fixed values of M and f v . Secondly, the 
larger value of R v reduces the outflow velocity-gradient such that 
the Sobolev optical depths of the line transitions become larger in 
the acceleration region. For our grid we considered R v -values of 3 
and 5 r min . 

As when making our comparison with the average spec- 
trum of Mrk 766 in Paper I, we constructed a multiplicative table 
("mtable") model for use with XSPEC (version 11, Arnaud 1996) 
from the ratio of our computed spectra to the primary power-law 
adopted in the radiative transfer simulations. We then attempted to 



fit the PG121 1+143 data using a model comprising of a power-law 
(with free normalization [iV] and power-law index [T]) combined 
with this mtable (which contains a total of 120 spectra from the 
grid of 3 x 2 x 20 in M, R v and /i). Given the strong but complex 
dependence of the spectral features to the line-of-sight (fi) we did 
not allow this parameter to be fit directly but stepped through the fi- 
values of the grid manually, freezing this parameter and then fitting 
the remaining four quantities (iV, Y, M and R v ). We did not allow 
for any component of systematic error in any of our fits as this is 
not expected to have a significant effect in the fit to the 2-10 keV 
region (see the discussion in Miller et al. 2009). Although other 
physical processes are expected to have a role in shaping the X-ray 
spectra of AGN (in particular, some degree of disk reflection and 
some warm absorber component seem all but inevitable), we did 
not include any additional component in this fit since our objective 
is to establish whether our disk wind models alone could be the 
dominant component that accounts for the major spectra features 
in the 2 - 10 keV band. 

For this model, the best fit was obtained for fi = 0.525 
(frozen), Y = 1.9, M = 0.68 M Q yr" 1 and R v = 3.2 r min 
and yielded xVd.o.f of 163/132. Although an imperfect fit, this 
suggests that the majority of the spectroscopic features are well ac- 
counted for by the model (a pure power-law fit to the same data 
yields x7 d -°- f of 389/134). 

To make a detailed comparison, we computed spectra for eigh- 
teen more models adopting parameters close to those favoured from 
the XSPEC fits and extending the range in R v (Y = [1.7, 1.9], 
R v = [3,5,8]r mln , M = [0.53,0.65,0.76] M yr" 1 ). For these 
simulations a larger number of MC quanta (~7x 10 6 ) were used 
to compute the spectra. The spectra obtained from each of these 
simulations are qualitatively similar although the amplitude of the 
spectra features varies slightly. We fit spectra from these eighteen 
models to the 2-10 keV spectrum in a similar manner to that de- 
scribed above - the best fitting spectrum obtained from this model 
grid (x 2 /d.o.f of 160/132) was for parameters close to that of the 
model with Y = 1.9, R v = 8r mm and M = 0.53 M Q yr _1 (for 
/i = 0.575); therefore we use the spectra from this model as the 
basis for our comparison (see below). 



5.3 Discussion 

Figure 5 compares the Fe Ka spectral region for the model with 
Rv = 8r mm , M = 0.53 M yr _1 and fi = 0.575 with the un- 
folded data of PG121 1+143 (Pounds & Reeves 2009). The com- 
parison is shown for the photon energy scale observed at the Earth, 
accounting for the source redshift of z = 0.0809 (Marziani et al. 
1996). 

The agreement between model and data is generally very good 
and shows that the disk wind paradigm readily produces the impor- 
tant features that characterize these data - blueshifted Fe Ka ab- 
sorption and broad red shifted Ka emission with a red-skewed wing 
- with strengths comparable to those observed. In agreement with 
the data, the model also predicts weak S XVI absorption around 
2.7 keV. There are, however, some clear discrepancies. For exam- 
ple, additional narrow components of emission may be required in 
the emission line profile (see discussion by Pounds & Reeves 2009) 
and the model seems to predict too little absorption at around 8 keV. 
However, none of these details suggest that the outflow paradigm 
is inappropriate for these data and may simply indicate physically 
distinct contributions to the spectrum. E.g. any very narrow emis- 
sion components required by the data are unlikely to be associated 
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Figure 5. Comparison of the Pounds et al. (2009) stacked data for 
PG121 1+143 (black crosses) and the wind model (black histogram) de- 
scribed in the text. The wind parameters are all identical to Model B (Ta- 
ble 2) except for V = 1.9, R v = 8r min and M = 0.53 M yr _1 . The 
model spectrum shown is for a viewing inclination of fi = 0.575. The red 
and blue histograms show the contributions to the model spectrum from 
direct and scattered/reprocessed radiation, respectively. The photon energy 
scale of the model spectrum has been adjusted to account for the redshift of 
PG121 1+143 (z = 0.0809). 

with the high velocity wind but can be more readily attributed to 
reflection by low-velocity material. 

Since the source is clearly complex, fitting the 2-10 keV 
spectrum to a wind model alone yields parameters that should only 
by regarded as indicative of those appropriate for PG121 1+143. 
To make a quantitative comparison of our models with the lower 
energy spectra for PG121 1+143 would requires significantly more 
complex fits of the data to be performed to address the strong soft 
excess and account for any physically distinct, lower-ionization 
contribution to absorption. However, at minimum, the highly- 
ionized wind should be able to imprint discreet blueshifted absorp- 
tion and broad emission features in the softer spectrum consistent 
with those that have been claimed from previous analyses of the 
PG121 1+143 RGS spectrum (e.g. Pounds et al. 2003; Pounds & 
Reeves 2007). 3 

To investigate this, we show in Figure 6 the 10 - 35A spec- 
tra (c.f. figures 6 - 9 of Pounds et al. 2003) for the same wind 
model used in Figure 5. The spectrum is shown for the line-of- 
sight (n — 0.575) that best fit the Fe K region of the stacked 
data of Pounds & Reeves (2009) and also for two higher inclina- 
tion angles {p = 0.475 and 0.375). The wavelength of blueshifted 
absorption lines for which identifications were claimed by Pounds 
et al. (2003) are indicated in the figure. Several of these - in par- 
ticular those corresponding to relatively high ionization states (e.g. 
Ne X, O VIII, N VII, C VI) - are clearly present in the model with 
strength and blueshift depending sensitively on the observer's in- 
clination. Fe XXIV imprints a clear additional absorption feature 
for jj, = 0.475 - although they did not discuss this in detail, we 
note that the fit in figure 9 of Pounds et al. (2003) does include an 

3 We note, however, that the exact combination of wind model and ob- 
server line-of-sight which fit the stacked EPIC pn data of Pounds & Reeves 
(2009) (which combines observations from 2001, 2004 and 2007) is not 
expected to be perfectly applicable for comparison with the RGS absorp- 
tion features discussed by Pounds et al. (2003) since it is known that the 
spectroscopic signatures of outflow are time-variable. 



Figure 6. The RGS region of the spectra of the same wind model shown in 
Figure 5. The black histogram shows the model spectrum for /i = 0.575 
(same inclination as used for the fit to the stacked EPIC pn data shown in 
Figure 5). The green and orange histograms show the spectra for fi = 0.475 
and 0.375, respectively. Each spectrum is shown as a ratio to the input pure 
power-law spectrum which had T = 1.9 and renormalised to 1.0 around 
30A. The identifications with red tick marks indicate the observed ener- 
gies of lines identified in the observations of PG121 1+143 by Pounds et al. 
(2003) that are distinct in the model spectra. The black tick marks indicate 
lines identified in the observations that do not have significant strength in 
any of the model spectra. Fe XXIV (blue tick mark) is responsible for the 
feature are 11 A in the fi = 0.475 spectrum which was not discussed by 
Pounds et al. (2003). 

absorption line around 1 1 A. The model also predicts significantly 
broadened O VIII Lyman a emission, qualitatively similar to that 
discussed by Pounds & Reeves (2007, 2009) in their more recent 
analyses of the PG121 1+143 RGS data. Thus our outflow models 
are able to simultaneously account for observable features in both 
the Fe K band and the softer-energy parts of the spectrum and there 
are good prospects that more detailed fitting of the complete set of 
spectral data will allow more robust constraints to be placed on the 
outflow model parameters. 

The less ionized species (Ne IX and O VII) are present in only 
very small quantities in the model and do not imprint lines in the 
spectra shown. The presence of these lines in the data suggests that 
a wider range of ionization states than in the model used here is 
required to explain the complete 2001 RGS spectrum. This may be 
found by exploring a greater range of our model parameter-space 
but may also require us to go beyond our assumption of a smooth, 
steady-state flow. A more realistic model incorporating flow inho- 
mogeneities will inevitably have a wider range of densities, and 
therefore ionization state. Such effects will be investigated in fu- 
ture studies. 

The spectral model used here is too simplistic to address the 
strong soft excess observed in PG121 1+143. However, we note that 
the wind model does clearly lead to an excess of emission below 
~ 1 keV suggesting that reprocessing by the fast flow may ac- 
count for a modest fraction of the excess soft-energy emission (c.f. 
King & Pounds 2003). Much of this contribution is due to heavily 
blended forests of emission lines (including those of L-shell Fe) but 
there are also moderately strong discreet features, such as the broad 
O VIII Lya emission line mentioned above. 

We defer more detailed studies including both a wider explo- 
ration of the model parameter space and the interplay of wind- 
formed spectral features with other components of the system to 
later work but the simple comparisons presented here already sup- 
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port the notion that the Fe K absorption and emission features 
could be predominantly formed in a fast outflow as in the picture 
presented by Pounds & Reeves (2009). Scaled to the black hole 
mass of PG1211+143 (4 x 10 7 M , Kaspi et al. 2000; see Sec- 
tion 2.8), the model with which we have compared has a mass- 
loss rate of M = 2.1 M yr _1 and a total covering fraction 4 of 
b = 0.7. These are comparable to the properties inferred by Pounds 
& Reeves (2009) and support their conclusion that the flow is very 
likely to be energetically important. 



6 CONCLUSIONS AND FUTURE WORK 

Outflows from the accretion disks around supermassive black holes 
are a promising explanation for the blueshifted absorption line fea- 
tures that have been identified in the X-ray spectra of AGN. They 
are also theoretically expected in high Eddington ratio sources for 
which line-driven winds have been modelled by e.g. (Proga & Kali- 
man 2004). AGN outflows may be massive enough to have signifi- 
cant implications for our understanding of accretion by supermas- 
sive black holes and likely have implications for the interpretation 
of a wide range of both absorption and emission features in the 
high energy spectra of AGN (Turner & Miller 2009). However, to 
properly interpret the observed spectroscopic features, modelling 
based on realistic theoretical spectra for outflow models is required 
in order to quantify the flow properties. Here, we have made a sig- 
nificant step towards this goal by extended the code described in 
Paper I to incorporate the physics of L- and M-shell ions and to 
compute both the ionization and thermal state of the outflow in de- 
tail. This allows us to explore a wider range of plausible outflow 
conditions including, in particular, less highly-ionized and more 
optically thick flows. 

An example calculation for one of our simply-parameterized 
wind models illustrates that, for mass-loss rates comparable to the 
Eddington accretion rate, a wide range of physical conditions are 
likely to be present in a fast outflow such that a very diverse range 
of spectral signatures are possible. As in Paper I, we find that con- 
siderable complexity can be introduced to the Fe Ka region where 
both narrow absorption and broad, red-skewed emission lines are 
predicted to form. However, such an outflow can also significantly 
affect the spectrum at softer energies - narrow, highly blueshifted 
absorption lines of lighter elements and lower ionization states of 
Fe are to be expected for many line-of-sight through an outflow. 
Emission features associated with these lines also form which can 
be broadened and skewed in a manner similar to the Fe Ka line, 
potentially allowing the wind model to account for broad Fe L fea- 
tures as have been reported in the spectrum of 1H0707-495 (see 
e.g. Fabian et al. 2009; Zoghbi et al. 2009). These features in the 
soft band are generally weaker that those in the Fe Ka region and 
their interpretation is more likely to be complicated owing to the 
wide range of physical components that may affect the softer bands 
(e.g. the soft excess and/or warm absorbers). 

As proof-of-concept, we have presented a simple compar- 
ison of our outflow model spectra with the well-known quasar 
PG121 1+143. We found that the most important features that have 
been identified in its hard X-ray spectra, namely a strong absorption 

4 Although our wind has a total covering fraction of b = 0.7, this value 
cannot be directly compared to that obtained from the much simpler model 
of Pounds & Reeves (2009) - there is significant diversity in the line-of- 
sight properties of our model and some of the lines-of-sight through the 
wind have only very low optical depth. 



line at ~ 7 keV and a broad emission feature peaking around 6 keV, 
can be simultaneously well-matched by our theoretical spectra. Al- 
though we have not fully explored the possible parameter-space of 
outflow models nor the interplay with other physically motivated 
phenomena (such as disk reflection), our results support the inter- 
pretation of Pounds & Reeves (2009) that the absorption and emis- 
sion components are likely physically related and form in a fast 
outflow that has a substantial covering fraction and a mass-loss rate 
comparable to the Eddington accretion rate. 

Our next step will be to attempt detailed model fits to one or 
more AGN using our radiative transfer code. This will allow us 
to more fully test the outflow paradigm and quantify the range of 
outflow properties which are suggested by observations, a critical 
step in establishing the role of highly-ionized outflows in the AGN 
phenomenon. 
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